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CONFIRMATION OF THE REMARKABLE COMPACTNESS OF MASSIVE QUIESCENT GALAXIES AT Z - 
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ABSTRACT 

Using deep near-infrared spectroscopy Kriek et al. (2006) found that ~ 45 % of massive galaxies at z ~ 2.3 
have evolved stellar populations and little or no ongoing star formation. Here we determine the sizes of these 
quiescent galaxies using deep, high-resolution images obtained with HST/NIC2 and laser guide star-assisted 
Keck/ AO. Considering that their median stellar mass is 1 .7 x 10 1 1 M Q the galaxies are remarkably small, with 
a median effective radius r e = 0.9 kpc. Galaxies of similar mass in the nearby Universe have sizes of w 5kpc 
and average stellar densities which are two orders of magnitude lower than the z ~ 2.3 galaxies. These results 
extend earlier work at z ~ 1.5 and confirm previous studies at z > 2 which lacked spectroscopic redshifts and 
imaging of sufficient resolution to resolve the galaxies. Our findings demonstrate that fully assembled early- 
type galaxies make up at most ~ 10 % of the population of ^-selected quiescent galaxies at z ~ 2.3, effectively 
ruling out simple monolithic models for their formation. The galaxies must evolve significantly after z ~ 2.3, 
through dry mergers or other processes, consistent with predictions from hierarchical models. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The sizes of massive galaxies at high redshift provide strong 
constraints on galaxy formation models. Models starting from 
ACDM initial conditions have predicted that the progenitors 
of today's early-type galaxies should be smaller by a fac- 
tor of a few at redshifts of 2-3 (e.g., Loeb & Peebles 2003; 
Robertson et al. 2006; Khochfar & Silk 2006; Naab et al. 
2007). This can be contrasted to expectations from simple 
"monolithic" models (e.g., Eggen, Lynden-Bell, & Sandage 
1962), in which galaxies undergo little or no structural evolu- 
tion since their initial assembly at very high redshift. 

Several recent studies have found evidence that massive 
galaxies at high redshift are indeed quite compact, particularly 
those with the lowest estimated star formation rates (Daddi 
et al. 2005; Trujillo et al. 2006, 2007; Zirm et al. 2007; Toft 
et al. 2007; Longhetti et al. 2007; Cimatti et al. 2008). The 
size evolution appears to be quite dramatic: Toft et al. (2007) 
find that the quiescent z ~ 2.5 galaxies in their sample have a 
factor of 30-40 higher surface density than red galaxies of the 
same mass in the Sloan Digital Sky Survey (SDSS). 

These results are still somewhat tentative, due to various 
systematic uncertainties. First, the galaxies were selected in 
relatively small fields, and therefore typically do not span the 
mass range of today's giant elliptical galaxies (M > 10 11 M ). 
Second, the imaging that has been used so far has either sam- 
pled the rest-frame ultraviolet (e.g., Cimatti et al. 2008), or 
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has not been of sufficient resolution to measure accurate rest- 
frame optical sizes of very compact objects. Third, some of 
the galaxies may have a central star burst or active nucleus 
which could skew their size measurements. Fourth, the stellar 
ages and masses of the galaxies have large uncertainties, as 
their broad-band photometry can typically be fitted by a wide 
range of models (see, e.g., Papovich, Dickinson, & Ferguson 
2001). Lastly, and perhaps most importantly, studies at z > 2 
are typically based on photometric redshifts, which are poorly 
calibrated for faint, red galaxies. 

Recently Kriek et al. (2006) used extensive near-infrared 
spectroscopy to secure the first sample of spectroscopically 
confirmed massive (M > 10 11 M©) galaxies with evolved stel- 
lar populations at z > 2. In this Letter, we present high qual- 
ity HST NIC2 and Keck images of these galaxies to measure 
accurate rest-frame optical sizes. The nine objects discussed 
here are the oldest objects in a complete sample of ^-selected 
massive galaxies at z ~ 2.3. If early-type galaxies were al- 
ready fully formed at this redshift we should see no struc- 
tural differences between these galaxies and galaxies of simi- 
lar mass at z = 0. 

2. SAMPLE, OBSERVATIONS, AND REDUCTION 

Kriek et al. (2006, 2008a) present a deep near-IR spec- 
troscopic survey of /if-selected galaxies at 2 < z < 2.7 with 
the Gemini Near-Infrared Spectrograph. The galaxies were 
originally identified in the deep, wide MUSYC near-IR imag- 
ing survey (Quadri et al. 2007). One of the key findings of 
Kriek et al. is that nearly half of the galaxies show no de- 
tected Ha emission and have spectra dominated by strong 
Balmer/4000 A breaks. Nine of these quiescent objects were 
presented in Kriek et al. (2006). These nine galaxies have a 
median redshift of 2.34, stellar masses > 10 11 M Q , and ages 
of 0.5 - 1 Gyr. They define a tight red sequence in rest-frame 
U-B (Kriek et al. 2008b). 

All nine quiescent galaxies were observed in the F160W 
filter using the NIC2 camera on HST, from June 2006 - June 
2007. Two orbits were used for each of the brightest two 
galaxies and three orbits for each of the remaining seven. 



1256-0 1256 - 1967 1256-142 



2 Sizes of Quiescent Galaxies at z ~ 2.3 

1030-1813 1030-2559 

Keck'- Keck K 



z = 2.56 



10 kpc 



r~\ — i — i — i — i — i — i — i — r 

4»V 



J , I i I i I , L 





~ i — 1 — i — 1 — i — I — i — 1 — r 



_i , i , i , i , L 



l — 1 — i — 1 — i — 1 — i — 1 — r 



J , I , I , I , L 



~i — 1 — i — 1 — i — 1 — i — 1 — r 



j , i , i , i_ 



0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 

ECDFS-5856 ECDFS- 11490 HDFS1-1849 HDFS2-2046 overage 



11490 HDFS1-1849 HDFS2-2046 overage 

• t • • 

z=2.34 z=2.31 z=2.24 <z>=2.34 




0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 0.3 0.4 0.5 0.6 0.7 

Fig. 1. — HST NIC2 images of the nine spectroscopically-confirmed quiescent z > 2 galaxies from Kriek et al. (2006). Each panel spans 
3 ."8 x 3 "8; North is up and East is to the left. The small panels below each galaxy show the best-fitting Sersic model (convolved with the PSF) 
and the residual after subtraction of the best-fitting model. The red ellipses are constructed from the best-fitting effective radii, axis ratios, and 
position angles. Note that the ellipses are significantly smaller than 10 kpc, which is the effective diameter of typical massive elliptical galaxies 
in the nearby Universe. Gemini GNIRS spectra from Kriek et al. (2006) are also shown. Insets show Keck LGS/AO images of three galaxies. 



Each orbit was split in two (dithered) exposures. The reduc- 
tion followed the procedures outlined in Bouwens & Illing- 
worth (2006), and R. Bouwens et al., in preparation. Before 
combining, the individual exposures were drizzled to a new 
grid with 0."0378 pixels to ensure that the PSF is well sam- 
pled. Images of the nine galaxies are shown in Fig.Q] 

Three of the galaxies (1030-1813, 1256-0, and 1256-1967) 
were also observed with Keck, using laser guide star assisted 
adaptive optics to correct for the atmosphere. The data were 
obtained on 2007 May 14 and 2008 January 13 using the 
NIRC2 wide field camera, which gives a pixel size of 0"04. 
The reduction followed standard procedures for near-IR imag- 
ing data. The Keck images are shown as insets in Fig. 1. They 
show the same qualitative features as the NICMOS data. 

3. FITTING 

Each galaxy was fitted with a Sersic (1968) radial surface 
brightness profile, using the 2D fitting code GALFIT (Peng 



et al. 2002). The Sersic n parameter allows for a large range 
of profile shapes, and provides a crude estimate of the bulge- 
to-disk ratio. For each galaxy a synthetic NIC2 PSF was cre- 
ated by generating subsampled PSFs with Tiny Tim 6.3 (Krist 
1995), shifting them to replicate the location of the galaxy on 
the individual exposures, binning these to the native NIC2 res- 
olution, and finally drizzling these "observations" to the grid 
of the galaxy images. The resulting fit parameters are listed 
in Table 1; ellipses corresponding to the best-fit parameters 
are indicated in red in Fig. Q] The (circularized) effective 
radii were transformed to kpc using Hq = 70 km s" 1 Mpc" 1 , 
fl m = 0.3, andO A = 0.7. 

Uncertainties in the structural parameters of faint galaxies 
are difficult to estimate, as they are usually dominated by sys- 
tematic effects. For each galaxy, we added the residual image 
of each of the other galaxies (excluding 1256-1967) in turn, 
repeated the fit, and determined the rms of the seven values 
obtained from these fits. The uncertainties listed in Table 1 
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Fig. 2. — Relations between size and (total) stellar mass (left panel) and between the average stellar density inside the effective radius and 
stellar mass (right panel). Large symbols with errorbars are the quiescent z ~ 2.3 galaxies. Small symbols are SDSS galaxies, with galaxies 
that are not on the red sequence in light grey. Broken lines indicate the expected location of galaxies with stellar velocity dispersions of 200, 
300, and 500 km s" 1 . The high redshift galaxies are much smaller and denser than SDSS galaxies of the same stellar mass. 



are 2x these rms values, to account for additional systematic 
uncertainties. These were assessed by changing the size of the 
fitting region, scrambling the subpixel positions of the galax- 
ies, and changing the drizzle grid. 



TABLE 1 

Structural Parameters 



ID 
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r d) 

' e 


± 


n 


± 


b/a 


± 


1030-1813 


2.56 


0.76 


0.06 


1.9 


0.5 


0.30 


0.03 


1030-2559 


2.39 


0.92 


0.18 


2.3 


0.6 


0.39 


0.04 


1256-0 


2.31 


0.78 


0.17 


3.2 


0.9 


0.71 


0.10 


1256-1967 


2.02 


1.89 


0.15 


3.4 


0.1 


0.75 


0.07 


1256-142 


2.37 


0.93 


0.04 


0.9 


0.3 


0.35 


0.04 


ECDFS-5856 


2.56 


1.42 


0.35 


4.5 


0.4 


0.83 


0.07 


ECDFS-11490 


2.34 


0.47 


0.03 


2.8 


0.8 


0.63 


0.07 


HDFS1-1849 


2.31 


2.38 


0.11 


0.5 


0.2 


0.29 


0.02 


HDFS2-2046 


2.24 


0.49 


0.02 


2.3 


0.8 


0.76 


0.08 



(1) Circularized effective radius in kpc. 



The Keck images offer an independent test of the reliability 
of the fit parameters. Fitting the Keck images with a range of 
stellar PSFs (including stars in the field-of-view) gives results 
that are consistent with the NIC2 fits within the listed uncer- 
tainties. As an example, for 1030-1813 we find r e = 0.73kpc, 
n = 1.6, and b/a = 0.32 from the Keck image. In the follow- 
ing, we will use the values derived from the higher S/N NIC2 
images; our conclusions do not change if we were to use the 
Keck results for 1030-1813, 1256-0, and 1256-1967. 

4. SIZES AND DENSITIES 

The most remarkable aspect of the z ~ 2.3 galaxies is their 
compactness. The circularized effective radii range from 0.5 
- 2.4 kpc, and the median is 0.9 kpc. To put this in con- 
text, this is smaller than many bulges of spiral galaxies (in- 
cluding the bulges of the Milky Way and M31, which have 
r e ss 2.5 kpc; van den Bergh 1999). In the left panel of Fig. [2] 
the sizes are compared to those of SDSS galaxies. The SDSS 
data were taken from the NYU Value Added Galaxy Catalog 



(Blanton et al. 2005) in a narrow redshift range, with various 
small corrections (see M. Franx et al., in preperation). Dark 
grey points are galaxies on the red sequence, here defined as 
u-g = 0.1 log(M) + (0.6 ±0.2). Stellar masses for the z ~ 2.3 
galaxies were taken from Kriek et al. (2008a) and corrected to 
a Kroupa (2001) IMF. The median mass of the z ~ 2.3 galax- 
ies is 1.7 x 10 11 Mq. The median r e of SDSS red sequence 
galaxies with masses 1.5—1.9 x 10 11 M Q is 5.0kpc, a factor 
of ~ 6 larger than the median size of the z ~ 2.3 galaxies. 

The combination of small sizes and high masses implies 
very high densities. The right panel of Fig. |2] shows the rela- 
tion between stellar density and stellar mass, with density de- 
fined as p = 0.5M /(|tt^) (i.e., the mean stellar density within 
the effective radius, assuming a constant stellar mass-to-light 
(M/L) ratio with radius). The median density of the z ~ 2.3 
galaxies is 3 x 10 10 M Q kpc -3 (with a considerable rms scatter 
of 0.7 dex), a factor of ~ 180 higher than the densities of local 
red sequence galaxies of the same mass. 

We note that it is difficult to determine the morphologies of 
the galaxies, as they are so small. Nevertheless, it is striking 
that several galaxies are quite elongated (see Fig.Q~|i. The most 
elongated galaxies are also the ones with the lowest n values 
(the correlation between n and b/a is formally significant at 
the > 99 % level 9 ), and a possible interpretation is that the 
light of a subset of the galaxies is dominated by very compact, 
massive disks (see §[5j . 

5. DISCUSSION 

We find that all (lOO^ %) of the quiescent, massive galax- 
ies at (z) = 2.3 spectroscopically identified by Kriek et al. 
(2006) are extremely compact, having a median effective ra- 
dius of only 0.9 kpc. This result extends previous work at 
Z ~ 1.5 (Trujillo et al. 2007; Longhetti et al. 2007; Cimatti 
et al. 2008), and confirms other studies at similar redshifts 
that were based on photometric redshifts and images of poorer 
quality (Zirm et al. 2007; Toft et al. 2007). Our study, to- 

9 There is no significant correlation between r e and n, or r e and b/a. 
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gether with the spectroscopy in Kriek et al. (2006) which 
demonstrates that the H-band light comes from evolved stars, 
shows that the small measured sizes of evolved high redshift 
galaxies are not caused by photometric redshift errors, AGN, 
dusty starbursts, or measurement errors. 

It is remarkable that all nine galaxies are so compact: even 
the largest galaxy in the sample (HDFS 1-1849) is signifi- 
cantly offset from the relations of red galaxies in the nearby 
Universe (see Fig. 2). We do not find any galaxy resembling 
a fully assembled elliptical or SO galaxy, which means that 
such objects make up less than ~ 10% of the population of 
quiescent galaxies at z ~ 2.3. This result effectively rules out 
simple monolithic models in which early-type galaxies are as- 
sembled at the same time as their stars, and is arguably the 
most direct evidence to date for an essentially hierarchical as- 
sembly history for massive galaxies. 

Our results are qualitatively consistent with recent hydro- 
dynamical simulations in a ACDM Universe, which predict 
that the stars in the central parts of massive elliptical galaxies 
formed rapidly at z ~ 5 out of a concentrated gas distribution 
(e.g., Naab et al. 2007). Observationally, it is tempting to link 
the compact z ~ 2.3 galaxies to the rapidly rotating, dense 
molecular gas reservoirs associated with some submm galax- 
ies 10 (Tacconi et al. 2008) and high redshift quasars (e.g., 
Riechers et al. 2006; Narayanan et al. 2008). The fact that 
several of the galaxies in our sample (and the two galaxies 
identified by Stockton et al. 2008) appear to be disks provides 
indirect support for this idea. 

Various mechanisms can play a role in bringing the z ~ 2.3 
galaxies onto the scaling relations followed by today's red 
galaxies. The total mass in galaxies on the red sequence in- 
creases by a factor of ~ 8 from z~2.3toz = (Kriek et al. 
2008b), and the galaxies that are added after z ~ 2.3 may have 
preferentially larger sizes than the galaxies that are already 
in place at that redshift. We note, however, that the galax- 
ies described here are quite extreme even when compared to 
the most compact 10 % of today's massive galaxies. Further- 
more, mass loss from stellar winds will decrease the stellar 
mass (and increase the radius) from z = 2.3 to z = 0. Finally, 
accretion of satellites (e.g., Naab et al. 2007) and dry mergers 
(e.g., van Dokkum 2005) may increase the sizes over time. 
Mergers will, to first order, move galaxies roughly parallel 
to the size - mass relation, but the effects on the sizes may 
be stronger in simulations with realistic boundary conditions 
(see Naab et al. 2007). Also, the smallest galaxies are likely 
to merge with larger galaxies, and even a merger between two 
small galaxies will (obviously) reduce their number. Each of 
these mechanisms could plausibly alter the size - mass rela- 
tion by a factor of 1.5-2, but not a factor of ~ 6. This means 

10 Note, however, that these submm galaxies would have to be at signifi- 
cantly higher redshift than the objects discussed in Tacconi et al. 



that some combination of effects is required to bring the com- 
pact z ~ 2.3 galaxies to the local relations — or that we have 
not yet identified the main mechanism. 

We also cannot exclude the possibility that we are overes- 
timating the densities of these galaxies, as there are several 
systematic uncertainties in the analysis. First, we may be un- 
derestimating the sizes of the objects if their morphologies 
are complex. In particular, an extended low surface bright- 
ness component could be "hiding" in the noise, in which case 
the measured effective radius would be larger if we had data 
of higher S/N. There is no evidence for this from fitting a Ser- 
sic profile to the stacked image of the nine galaxies (see Fig. 
[T]), but we do note that the galaxies are on average ~ 0.2 mag 
brighter in ground-based images (which include any extended 
emission) than in the HST images. A more subtle effect is an 
age-induced radial gradient in the M/L ratio of the galaxies, 
such that the central parts are younger than the outer parts. If 
this is the case the mass-weighted r e could be significantly 
larger than the luminosity-weighted r e . Such a spatial ar- 
rangement arises in several different classes of models for the 
formation of massive ellipticals (Robertson et al. 2006; Naab 
et al. 2007, R Hopkins et al., in preparation). Lastly, the IMF 
might be weighted toward high mass stars in these systems. 
Several studies have argued that the IMF in the progenitors of 
elliptical galaxies may have been deficient in low mass stars, 
which would impact their inferred masses and densities (e.g., 
Larson 2005, van Dokkum 2008). The precise effect depends 
on the assumed form of the IMF, and can be a factor of ~ 2 
for plausible parameters (see van Dokkum 2008). 

There are several ways to improve upon the results pre- 
sented here. Radial gradients in the M/L ratio can, in prin- 
ciple, be measured from high quality imaging in multiple 
passbands. Very deep imaging can provide better constraints 
on the surface brightness profiles at large radii. Ultimately, 
kinematics should provide the most direct information on the 
nature of these compact objects (see also Toft et al. 2007, 
Cimatti et al. 2008). As shown by the broken lines in Fig. 
2 (derived from the relation between M, a, and r e given in 
van Dokkum & Stanford 2003) the implied velocities are very 
high at 300 - 500 km s" 1 . Finally, it will be interesting to build 
up spectroscopic samples that are fainter in K, as this will al- 
low us to extend the analysis to lower masses and to determine 
the sizes of quiescent massive galaxies with higher M/L ratios 
than the objects studied here. 
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